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Abstract 

We compute the leading order (in ag) perturbative QCD and power (1/m^) corrections to 
the hadronic invariant mass and hadron energy spectra in the decay B — s- Xsi'^£~ in standard 
model. The computations are carried out using the heavy quark expansion technique (HQET) 
and a perturbative-QCD improved Fermi motion (FM) model which takes into account B-meson 
wave-function effects. The corrections in the hadron energy (Eh) spectrum are found to be small 
over a good part of this spectrum in both methods. However, the expansion in 1/mb in HQET 
fails near the lower kinematic end-point and at the cc threshold. The hadronic invariant mass 
(Sh) spectrum is calculable only over a limited range Sh > Ams in the heavy quark expansion, 
where A ~ ms — mf,. We also present results for the first two hadronic moments {S^) and {E'^), 
n — 1,2, working out their sensitivity on the HQET and FM model parameters. For equivalent 
values of these parameters, the moments in these methods are remarkably close to each other. The 
constraints following from assumed values of {S^) on the HQET parameters Ai and A are worked 
out. Data from the forthcoming B facilities could be used to measure the short-distance contribution 
in B Xs£+i!"and constrain the HQET parameters Ai and A. This could be combined with 
complementary constraints from the decay B ~* Xivi to determine these parameters precisely. We 
also study the effect of the experimental cuts, used recently by the CLEO collaboration in searching 
for the decay B Xsi^i~, on the branching ratios, hadron spectra and hadronic invariant mass 
moments using the FM model. 
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1 Introduction 



The semileptonic inclusive decays B Xgi^i^ , where = e^,fi^,T^, offer, together with the 
radiative electromagnetic penguin decay B Xg + 7, presently the most popular testing grounds for 
the standard model (SM) in the flavor sector. This is reflected by the impressive experimental and 
theoretical activity in this field, reviewed recently in and [^, respectively. We shall concentrate 
here on the decay B — > X.si'^i^ for which the first theoretical calculations were reported a decade 
ago [§-|5|, emphasizing the sensitivity of the dilepton mass spectrum and decay rate to the top quark 
mass in the short-distance contribution. With the discovery of the top quark and a fairly accurate 
measurement of its mass theoretical emphasis has changed from predicting the top quark mass 
using this decay to using its measured value as input and making theoretically accurate predictions 
for the decay rates and spectra. This will help confront the predictions in the SM with experiment 
more precisely and will allow to search for new phenomena, such as supersymmetry [f7|-pT|]. 

Since these early papers, considerable theoretical work has been done on the decay B Xgi^i^ in 
the context of the standard model. This includes, among other aspects, the calculation of the complete 
leading order perturbative corrections in the QCD coupling constant Og to the dilepton invariant 
mass spectrum forward-backward (FB) asymmetry of the leptons |14,1^, and, additionally, 

leading order power corrections in l/m^ to the decay rate, dilepton invariant mass spectrum and FB 
asymmetry [^] , using the heavy quark expansion technique (HQET) We recah that the 

corrections to the dilepton spectrum and decay rate in B ^ Xsi~^£~weice calculated in ref. |jl8| but 
their results were at variance with the ones derived later in ref. The power corrected dilepton 

mass spectrum and FB asymmetry have been rederived for the massless s-quark case recently [|ig|], 
confirming the results in ref. |15|. Corrections of order to the dilepton mass spectrum away from 
the {J/tpji/^' , ...)-icesonant regions have also been worked out [2^,21|, making use of earlier work on 
similar power corrections in the decay rate for B Xg +7 p2,23|. The l/m^ power corrections to the 
left-right asymmetry |24,25| have been presented in |jl^ correcting an earlier calculation of the same 
[p5|]. Likewise, the longitudinal polarization of the lepton, Pl, in i? — > XsT~^t~ at the partonic level 
has been worked out ^6[; the other two orthogonal polarization components Pt (the component in 
the decay plane) and Pl (the component normal to the decay plane) were subsequently worked out in 
ref. 1^^. As an alternative to HQET, S-meson wave-function effects in the decay B — > Xsi~^(~ have 
also been studied for the dilepton invariant mass spectrum and FB asymmetry jl^, using the Fermi 
motion (FM) model [p8|. Some of the cited works have also addressed the long-distance aspect of the 
decay B — > X^^'''^" having to do with the resonant structure of the dilepton invariant mass spectrum. 
We shall leave out the J/i/j, -0', ...-resonant contributions in this paper and will present a detailed 



phenomenological study including them elsewhere |29]. 
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This theoretical work, despite some uncertainties associated with the LD-part, will undoubtedly 
contribute significantly to a meaningful comparison of the SM and experiment in the decay B — > 
Xgi'^i". Still, concerning the SD-contribution, some aspects of this decay remain to be studied 
theoretically. In the context of experimental searches for B — > Xsi'^i~ , it has been emphasized (see, 
for example, the CLEO paper |^^) that theoretical estimates of the hadronic invariant mass and hadron 
energy spectra in this decay will greatly help in providing improved control of the signal and they will 
also be needed to correct for the experimental acceptance. In addition to their experimental utility, 
hadron spectra in heavy hadron decays are also of considerable theoretical interest in their own right, 
as reflected by similar studies done for the charged current induced semileptonic decays B Xcii^i 



and B Xu£i/£ [31-35], where the main emphasis has been on testing HQET and/or in determining 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements Vet and Vui,. The hadronic invariant mass 
spectra in 6 ^ si~^i~ and b ui~U£ decays have striking similarities and differences. For example, 
both of these processes have at the parton level a delta function behavior dT/dsQ oc 6{sQ—mg), q = u,s, 
where sq is the hadronic invariant mass at the parton level. Thus, the entire invariant mass spectrum 
away from sq = is generated perturbatively (by gluon bremsstrahlung) and through the S-hadron 
non-perturbative effects. Hence, measurements of these spectra would lead to direct information on 
the QCD dynamics and to a better determination of the non-perturbative parameters. There are also 
obvious differences in these decays, namely the decay B — > Xuii^e is intrinsically a lot simpler due to 
the absence of the resonating cc contributions, which one must include to get the inclusive spectra 
in i? — > Xsi~^£~, or else use data in restricted phase space where the cc-resonant contributions are 
subleading. 

Having stated the motivations, we study hadron spectra in the decay B Xsi'^i~ in this pa- 
per. We first compute the leading order (in Os) perturbative QCD and power (1/m^) corrections 
to the hadronic invariant mass and hadron energy spectra at the parton level. In addition to the 
bremsstrahlung contribution 6 — > (s -|- g)i^i^ , there are important non-perturbative effects even in 
0(0°) that come from the relations between the b quark mass and the B meson mass. In HQET, this 
takes the form niB = nib -|- A — (Ai -|- 3X2) /2mb + where A, Ai and A2 are the HQET parameters 
[p^-p!8[. Keeping, for the sake of simplicity just the A term, the hadronic invariant mass Sh is related 
to So and the partonic energy Eq by Sh = •50 + 2A£'o + A^. This gives rise to a non-trivial spectrum 
in the entire region A^ < Sh < M"^. Including both the 0{l/rn^) and 0{as) terms generates hadron 
energy and hadronic invariant mas spectrum with terms of 0{A/mB), 0{asA/mB), 0{\i/m\) and 
0{\2/m?B)- '^^^ power- and perturbatively corrected hadron spectra up to and including these terms 
are presented here. The 1/m^ corrections in the hadron energy spectrum are found to be small over a 
good part of this spectrum. However, the expansion in 1/mf, fails near the lower end-point and near 
the cc threshold. The hadronic invariant mass spectrum is reliably calculable over a limited region 
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only, namely for Sh > Ams. Hadronic moments (5"]^) and {E'^), on the other hand, are calculable in 
HQET and we have summarized the results for the first two moments n = 1, 2 in a letter |^6[, based 
on this study. The hadronic invariant mass moments are sensitive to the HQET parameters A and 
Ai. This provides potentially an independent determination of these quantities. We think that the 
hadron spectra in 5 ^ Xs^^£~and B X^ii^i can be related to each other over limited phase space 
and this could help in vastly improving the present precision on Vub [0 and the parameters Ai and A 

mm- 

In view of the continued phenomenological interest in the FM model |28|, motivated in part 



by its close resemblance to the HQET framework ||39|,17], we also compute the hadron spectra in 
B — > Xs£'^£~ in this model, taking into account the O^Og) perturbative QCD corrections. The FM 
model is characterized by two parameters which are usually taken as pp, the Gaussian width of 
the 6-quark momentum, and rriq, the spectator quark mass in the B hadron; the 6-quark mass is a 
momentum-dependent quantity (see section 6 for details). The matrix element of the kinetic energy 
operator, Ai and the binding energy A can be calculated in terms of the FM model parameters. The 
difference between the effective 6-quark mass, which is a derived quantity in the FM model, and the 
5-meson mass can also be expressed via an HQET-type relation, nis = m^^+A — Ai/2m^^. However, 
there is no analog of A2 in the FM model. Having defined the equivalence between the FM model 
and HQET parameters, we shall use A and Ai to also characterize the FM model parameters. The 
dependence of the hadron spectra in the FM model in the decay B XsH.'^l'^ on the parameters A 
and Ai is studied in this paper. We find that the hadron energy spectrum in i? — > Xs^^£~in the FM 
model is stable against variations of the model parameters. The hadron energy spectra in the FM 
model and HQET are also found to be close to each other in regions where HQET holds. This feature 



was also noticed in the context of the decay B X^ii^e in ref. |33]. The hadronic invariant mass 
spectrum depends sensitively on the parameters of the FM model - a behavior which has again its 
parallel in studies related to the decay B Xuii^e as well as in HQET. Hadronic moments (5]^) 
and {E^) are computed in the FM model and are found to be remarkably close to their counterparts 
calculated in HQET for equivalent values of the parameters A and Ai. The picture that emerges from 
these comparisons is that the spectra and moments in the two approaches are rather similar, though 
not identical. We also study the effects of the CLEO experimental cuts on the branching ratios, hadron 
spectra and hadronic moments in B ^ Xgi'^l" in the FM model and the results are presented here. 
These can be compared with data when they become available. 

This paper is organized as follows. In section 2, we define the kinematics of the process B 
Xs£~^i~ and introduce the quantities of dynamical interest in the framework of an effective Hamil- 
tonian. Leading order (in Os) perturbative corrections to the hadron energy and hadronic invariant 
mass spectra at the parton level are derived in section 3, where we also present the Sudakov-improved 
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spectrum dB/dsQ. Using the HQET relation between ms and m^, we calculate the corrected hadronic 
invariant mass spectrum dB/dSu- In section 4, we present the leading power corrections (in 
for the Dalitz distribution d^B / dx^dsQ (here xq and sq are the scaled partonic energy and hadronic 
invariant mass, respectively) and derive analytic expressions for the hadron energy spectrum dB/dxQ 
and the resulting spectrum is compared with the one in the parton model. In section 5, we calculate 
the moments in the hadron energy and hadronic invariant mass in HQET and give the results for 
(Sh), {Sjj), {Eh) and {Ej^) in terms of the corresponding moments in the partonic variables. Sec- 
tion 6 describes the wave-function effects in the FM model |^8[ in the hadron energy and hadronic 
invariant mass spectra. We also give here numerical estimates of the hadronic moments in HQET and 
the FM model. In section 7, we study the effects of the experimental cuts used in the CLEO analysis 
of B ^ Xgi^i^on the hadron spectra and hadronic moments using the FM model. Estimates of the 
branching ratios B{B Xgi'^i~) for i = fi,e are also presented here, together with estimates of the 
survival probability for the CLEO cuts, using the FM model. Section 8 contains a summary of our 
work and some concluding remarks. Definitions of various auxiliary functions and lengthy expressions 
appearing in the derivation of our results, including the partonic moments {xq), {{sq — rhs)'^) and 
(xo(so — w-s)) for n = 1, 2 are relegated to the Appendices A - D. 

2 The Decay B Xst^l^ in the Effective Hamiltonian Approach 
2.1 Kinematics 

We start with the definition of the kinematics of the decay at the parton level, 

h{Pb) ^ s{ps){+g{pg)) + (.^{p+)+r{p-) , (1) 

where g denotes a gluon from the 0{as) correction (see Fig. 1). The corresponding kinematics at the 
hadron level can be written as: 

B{pb) ^ Xs{ph) + ^+(P+) + r (p_) . (2) 

We define the momentum transfer to the lepton pair and the invariant mass of the dilepton system, 
respectively, as 

q = p++p_ , (3) 
s = q^. (4) 

The dimensionless variables with a hat are related to the dimensionful variables by the scale rrib, the 
6-quark mass, e.g., s = rhg = — etc.. Further, we define a 4- vector t>, which denotes the velocity 



4 



of both the 6-quark and the S-meson, pb = rribV and pB = tubv. We shall also need the variable u 
and the scaled variable n = -t-, defined as: 

D 

u = -{pb - p+f + {pb - P-f , (5) 
u = 2v ■ — p_) . (6) 

The hadronic invariant mass is denoted by Sh = P% and Eh denotes the hadron energy in the final 
state. The corresponding quantities at parton level are the invariant mass sq and the scaled parton 
energy = In parton model without gluon bremsstrahlung, this simplifies to sq = Tn^ and 

xq becomes directly related to the dilepton invariant mass xq = 1/2(1 — s + "fhl). From momentum 
conservation the following equalities hold in the 6-quark, equivalently i?-meson, rest frame {v = 
(1,0,0,0)): 

xo = l- v- q, so = \-2v-q + s, (7) 
Eh = niB — V ■ q , Sh = — 2mBV ■ q + s . (8) 

The relations between the kinematic variables of the parton model and the hadronic states , using the 
HQET mass relation, can be written as 

- A1 + 3A2 / V Ai + 3A2\ 

Eh = A +{mB-A + — Uq + . . . , 

ZniB V ZniB / 

Sh = m1 + P + (m| - 2AmB + + Ai + 3A2) (sq - m^) 

+ {2AmB - 2A^ - Ai - 3A2)a;o + . . . , (9) 

where the ellipses denote terms higher order in l/nib. 
2.2 Matrix element for the decay B Xsi^£~ 

The effective Hamiltonian obtained by integrating out the top quark and the bosons is given as 

6 

Ci{pL)Oi + Cj{^i)-^Saa^,u{mbR + msL)baF'"' 



.1=1 



,(10) 



where L and R denote chiral projections, L{R) = 1/2(1 =p 75), Vij are the CKM matrix elements and 
the CKM unitarity has been used in factoring out the product V*gVtb- The operator basis is taken 



from 1 15], where also the Four-Fermi operators Oi, . . . ,0q and the chromo-magnetic operator Og can 
be seen. Note that Og does not contribute to the decay B Xgi^i^ in the approximation which 
we use here. The Ci{^) are the Wilson coefficients, which depend, in general, on the renormalization 
scale /i, except for Ciq. 



The matrix element for the decay B — > Xs£~^i can be factorized into a leptonic and a hadronic 
part as 



with 



M{B ^Xsi+t 



^ 7/. {cf{s) T Cio + 2C7«ff 4\+2ms cf 7m 4^ 



6. 



(11) 

(12) 
(13) 



rr 

The effective Wilson coefficient Cg (5) receives contributions from various pieces. The resonant cc 

rr 

states also contribute to Cg (5); hence the contribution given below is just the perturbative part: 



cf{s) = C,rj{s) + Y{s) 



(14) 



Here rj{s) and Y[s) represent the 0{as) correction [40| and the one loop matrix element of the Four- 
Fermi operators p^jl^ , respectively. While Cg is a renormalization scheme-dependent quantity, this 
dependence cancels out with the corresponding one in the function Y{s) (the value of ^, see below). 

rr 

To be self-contained, we list the two functions in Cg (s): 



Y{s) = g{m^,s){'iCi + C2 + ^C^ + Ci + ?,C^ + CQ) 

-iff(l, s) (4 C3 + 4 Q + 3 C5 + Ce) - ^5(0, s) (C3 + 3 

+1 (3C3 + C4 + 3C5 + Q) - (3Ci + C2 - Cs - 3C4) , 

rj{s) = IH uj(s) , 



(NDR), 
-1 (HV), 



(15) 
(16) 

(17) 



5(0, s) = 
where y = 4z^/s, and 
a;(s) 



8 rrib 

-y^^^T^ 

G(l - y)(ln ^ ^ ~ ^ - ivr) + @(y - 1)2 arctan ^ ^ 

1 - VI -y Vy - 1. 

8 8 mfc 4 4. 

ml — ) m s H — zvr , 

27 9 ^ // ^ 9 9 



2 2 4 

TT 1^^2 S - 

9 3 ^ ^ 

2s(l-Fs)(l-2s) 
3(l-s)2(l-K2s) 



2 5 -|- 4s 

- In s ln(l - s) - ln(l - s) 

3 3(1 -|- 2s) 



In s -|- 



5 + 9s- 6s2 
6(l-s)(l-F2s) 



(18) 
(19) 



(20) 
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Above, (NDR) and (HV) correspond to the naive dimensional regularization and the 't Hooft-Veltman 
schemes, respectively. The one gluon correction to O9 with respect to xq will be presented below in 



eq. (26). The Wilson coefficients in leading logarithmic approximation can be seen in [[T^ ]. 

With the help of the above expressions, the differential decay width becomes on using p± = 
iE±,p±), 

dr = \VZV,,f /'f+ /'f- (W\^ L^"' + L^"-") , (21) 

2mB 27r2 ' *^ (27r)32^+ (27r)32^_ V i^^y j , \ j 

where W^^^ and Lj^jJ^ are the hadronic and leptonic tensors, respectively. The hadronic tensor w/j^J^ 
is related to the discontinuity in the forward scattering amplitude, denoted by T^J, , through the 
relation Wf^i, = 21mTfj_y. Transforming the integration variables to s, u and v ■ q, one can express the 
Dalitz distribution in B ^ X^i^i~as: 



auasa{v ■ q) 2mB 2 7r"= 256 vr* ^ 

with 



(22) 



r^/V = i/d^ye-*^-'^(5|T{ri^/^(y),r2^/^(0)}|5) , (23) 



L^/^'^" = 2 



L/r\ _ -p L/R _ j^L/R 



where Fi^ = = T^ , and is given in eq. (13). Using Lorentz decomposition, the tensor 



can be expanded in terms of three structure functions Tj, 

T^u = -Ti g^u + T2 Vy + Ts ie^^ap , (25) 

where the structure functions which do not contribute to the amplitude in the limit of massless leptons 
have been neglected. The problem remaining is now to determine the Tj, to which we shall return in 
section Q 

3 Perturbative QCD Corrections in 0{as) in the Decay B Xsi~^i~ 

In this section the 0{as) corrections to the hadron spectra are investigated. Only Og is subject to 
corrections and the renormalization group improved perturbation series for Cg is 0{l/as) + 0(1) + 
0{as) + . . ., due to the large logarithm in Cg represented by 0{l/as) [0]. The Feynman diagrams, 
which contribute to the matrix element of O9 in 0(qs), corresponding to the virtual one-gluon and 
bremsstrahlung corrections, are shown in Fig. |l|. The effect of a finite s-quark mass on the 0{as) 
correction function is found to be very small. After showing this, we have neglected the s-quark mass 
in the numerical calculations of the 0{as) terms. 
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Parameter 


Value 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


sm Ow 


0.2325 


rris 


0.2 (GeV) 


rric 


1.4 (GeV) 


rrib 


4.8 (GeV) 


nit 


lib ± 5 (GeV j 










129 


as{mz) 


0.117 ±0.005 


Bsi 


(10.4 ± 0.4) % 



Table 1: Default values of the input parameters and errors used in the numerical calculations. 



Ci 


C2 




C4 








Cg 


Cio 


C(o) 


-0.240 


+1.103 


+0.011 


-0.025 


+0.007 


-0.030 


-0.311 


+4.153 


-4.546 


+0.381 



Table 2: Values of the Wilson coefficients used in the numerical calculations correspondinq to the 
centra, „lue, oj tke .ara^eter, ,.,en ,n TaMe Here, cf . Cr - ft/3 - ft. a„, for ft ™ 
the NDR scheme. 

3.1 Hadron energy spectrum 

The explicit order as correction to Og can be obtained by using the existing results in the literature as 
follows: The vector current Og can be decomposed as V = {V — A)/2 + (V + A)/2. We recall that the 
(V — A) and {V + A) currents yield the same hadron energy spectrum [^] and there is no interference 
term present in this spectrum for massless leptons. So, the correction for the vector current case in 



B Xsi^i can be taken from the corresponding result for the charged (V — A) case [28,40|, yielding 

(26) 



Cfixo) = Cgpixo) + Y{xo) 




Figure 1: Feynman diagrams contributing to the explicit order corrections of the operator Og. 
Curly lines denote a gluon. Wave function corrections are not shown. 
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with 



p{x) 
(t{x) 



TT 



(27) 
(28) 



(3x — — 2m^ + Sm^x) S^x^ — ' 
where i^(xo) = i^(s) with s = 1 — 2xo+m^. The expression for Gi{x) with 7^ has been calculated 
in [40|. The effect of a finite rUg is negligible in Gi(x), as can be seen in Fig. |2|, where this function 
is plotted both with a finite s-quark mass, rus = 0.2 GeV, and for the massless case, = 0. A 
numerical difference occurs at the lowest order end point x™"^ = 1/2(1 + m^) (for mi = 0), where 
the function develops a singularity from above (xo > x^""^) and the position of which depends on the 
value of rris- The function Gi(x) for a massless s-quark is given and discussed below [40|. 



Gi(x) 



Gi(x) 



x'^{ — {lQx^ - 84x^ + 585x2 _ -^ggg^ ^ -^215) + (8x - 9) ln(2x) 
90 



+ 2(4x - 3) 



TT 



+ Lz2(l-2x) 



} for < X < 1/2 , 



Y^(l - x)(32x^ - 136x^ + 1034x3 - 2946x2 + 1899x + 312) 
^ ln(2x - l)(64x3 - 48x2 - 24x - 5) 



+ x^(3-4x) 
where Li2{z) is the dilogarithmic function 



^ - iLi2{^) + ln2(2x - 1) - 21n2(2x) 
3 2x 



for 1/2 < X < 1 



(29) 



0.5 



-0.5 




Figure 2: The function Gi{x) is shown for rUs = 0.2 GeV (solid line) and for the massless case 
corresponding to eq. (dashed line). 



The 0{as) correction has a double logarithmic (integrable) singularity for xq 1/2 from above 
(xo > 1/2). Further, the value of the order corrected Wilson coefficient Cg^(xo) is reduced 



compared to its value with = 0, therefore also the hadron energy spectrum is reduced after including 
the explicit order QCD correction for < xq < 1/2. Note that the hadron energy spectrum for 
B Xs^"''^" receives contributions for 1 > x > 1/2 only from the order bremsstrahlung corrections. 

3.2 Hadronic invariant mass spectrum 
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Figure 3: The differential branching ratio ^ parton model is shown in the 0{as) 

bremsstrahlung region. The dotted (solid) line corresponds to eq. (eq. (jg^j ). The vertical line 

denotes the one particle pole from b s£~^£~ . We do not show the full spectra in the range < sq < ni^ 
as they tend to zero for larger values of sq . 

We have calculated the order ag perturbative QCD correction for the hadronic invariant mass in 
the range < sq < 1. Since the decay 6 —> s + £~ contributes in the parton model only at 
So = rhi, only the bremsstrahlung graphs 6— > s + g + + l~ contribute in this range. This makes 
the calculation much simpler than in the full sq range including virtual gluon diagrams. We find 

^ = ^^0— j-{ ^''°~'^\ 93 - 41so - 95s2 + bbsl) + ^ lnso(-3 - 5so + 9% - 2s^)}C| . (30) 
dso 6 t: sq z( y 

Our result for the spectrum in 5 ^ X^^^^^is in agreement with the corresponding result for the 
(y — A) current obtained for the decay B — > Xqlv^ in the ruq = limit in ref. (their eq. (3.8)), 
once one takes into account the difference in the normalizations. We display the hadronic invariant 
mass distribution in Fig. |3| as a function of sq (with sq = tu^sq), where we also show the Sudakov 
improved spectrum, obtained from the 0{as) spectrum in which the double logarithms have been 
resummed. For the decay B X^iui, this has been derived in ref. |33], where all further details 
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Figure 4: The differential branching ratio ^^^^"j^^ ^ ^ ^/^e hadronic invariant mass, Sh, shown 
for different values of rub in the range where only hremsstrahlung diagrams contribute. We do not 
show the result in the full kinematic range as the spectra tend monotonically to zero for larger values 
of Sh < m|. 



can be seen. We confirm eq. (17) of ref. |Q for the Sudakov exponentiated double differential decay 
rate g^j^ and use it after changing the normalization Tq Bo^Cg for the decay B Xsi'^i^ ■ The 
constant Bq is given later. Defining the kinematic variables (x, y) as 



2 2 

X ruf^ , 
1 



v-q = {x + ^{l - xfy)mb . 



the Sudakov-improved Dalitz distribution is given by 
d'^B 



dxdy 



{B^Xst+n = -Bo^x{l-x^f{l + 2x')exp[-'^ln\l-y) 



Aos ln(l - y) 
37T (l-y) 



2a, 



l-^iGix)+Hiy)) 



2as dH 
Sir dy 



(y) cl , 



where [|3^ 



(31) 



(32) 



Gix) 



[8x^(1 -x^ - 2x^) \nx + 2(1 - x^)^(5 + Ax^) ln(l - x^) - (1 - x^)(5 + - Qx'^ 



2(l-x2)2(l + 2x2) 



W + 2Li2{x^) - 2Li2(l - x^) 
2 - z{l - x) + K 



(33) 



H{y) = / dz(- In. 

Jo ^ 1 — 2; 



[1 — x){3 + X + xz — z) 

(1 + X)2 



ln(l - z) - 2 In 
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2 - z(l - x) + K 



7(l + x)(l + 2x^^^^_^^^3-2x^)l). (34) 



2(1 +x)2(l + 2x2) L 1 
The quantity k in eq. (|3^ is defined as k = \/z'^{l — x)^ + Axz. 

To get the hadronic invariant mass spectrum for a b quark decaying at rest we change variables 
from {x,y) to (q'^jSq) followed by an integration over q'^, 

dB _ r("^b~V^f (Pb 1 
dso J4mf dxdy 2m^x(l - x)^ 

The most significant effect of the bound state is the difference between uib and nib, which is 
dominated by A. Neglecting Ai, A2, i.e., using A = uib — mb, the spectrum is obtained along 

dB 

dso ' 

rriB —A 



the lines as given above for after changing variables from (x,y) to (q^jSH) and performing an 
integration over q'^. It is valid in the region 771^—^^—3^^ < Sh < w-^ (or < Sh < "i-^, 



dB__ ri^B-ysJif 1 
d5_ff jAmf dxdy 2mfmBx(l - x) 



neglecting rris) which excludes the zeroth order and virtual gluon kinematics (sq = ^1)-, yielding 

dq^^ L . (36) 

The hadronic invariant mass spectrum thus found depends rather sensitively on rrih (or equivalently 
A), as can be seen from Fig. ^. An analogous analysis for the charged current semileptonic B decays 
B Xuii'e has been performed in ref. [B4|, with similar conclusions. 



4 Power Corrections in the Decay B Xsi^t 



The hadronic tensor in eq. ( |25D can be expanded in inverse powers of ruh with the help of the HQET 
techniques. The leading term in this expansion, i.e., 0{m^) reproduces the parton model result. In 
HQET, the next to leading power corrections are parameterized in terms of the matrix elements of 
the kinetic energy and the magnetic moment operators Ai and A2, respectively. The B — B* mass 
difference yields the value A2 = 0.12 GeV^. In all numerical estimates we shall use this value of A2 and, 
unless otherwise stated, we take the value for Ai extracted from an analysis of data on semileptonic 
B-decays {B Xii^i), yielding Ai = -0.20 GeV^ with a corresponding value A = 0.39 GeV For 



a review on the dispersion in the present values of these non-perturbative parameters, see [38| 



The contributions of the power corrections to the structure functions Tj can be decomposed into 
! sum of various terms, denoted by T- , which 
evaluation of the time-ordered product in eq. (p3|): 



(i) 

the sum of various terms, denoted by T- , which can be traced back to well defined pieces in the 



Tiiv.q,s)= J2 T^^'\^-Q,s). (37) 

j=0,l,2,s,g,5 

The expressions for Tp\v.q, s), i = 1,2,3 calculated up to 0{mB/ml) are given in [|l^]. After 
contracting the hadronic and leptonic tensors, one finds 
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With the help of the kinematic identities given in eq. (|^) , we can make the dependence on xq and sq 
exphcit, 



r^/^^, L^/^''" = rrib^ 1 2(1 - 2xo + so)ri^/^ + - ^u^ - sq Ts^/^ ^ (1 - 2xo + so)n T- 



(39) 



and with this we are able to derive the double differential power corrected spectrum ^^^^^^ for B 
XJ+i~. Integrating eq. (PD over u first, where the variable u is bounded by 



- 2^x'o -so<u< +2^x'o - so , (40) 
we arrive at the following expression 

- ^BolmJxl - so I (1 - 2X0 + so)Ti [sq, xq) + ^^^T2{so, xq) 1 + ©(A^a,) , (41) 



dxQ dso 
where 



Ti(so,xo) = ^|(^8xo-4(^ + A2)^ (|C|ff(l)|2 + |Cio|2 



+ ^32(-2m^ - 2so - 4m^so + xq + 5m^xo + soXq + m^soXo) + 16(^ + A2) 

|(^eff|2 

X (-5 - llm^ + 5so - "^po + lOxo + 22m^xo - IOxq - lOm^xg) > ' 



(so - 2x0 + 1)2 



+ 



-32 



.(^2 + _ ^„ _ ^2^q) _ 48(^ + A2) ) Reicfis)) Cf] 



so-2xo + r ' ' ' "3 7 

+ ^ I (^(-2^0 - 3x0 + 5x2) ^ 8A2(-2so + xo + Sxg)^ (|C|^(s)|2 + jCio'^^ 

(32^]^ ,0 9999 9 9 

— ^ — (Sm^ + 12so + ISm^so — 2so — 2m^So — 3xo — 217713X0 — 13soXo — lOm^soXo 

— 3xq + Om^xQ + SsqXq + Stti^sqXq + 4xo + 4mgXQ) 

+ 32A2(— 2m2 — 2m1so — 2sq — ^in^sl + xo — rn^xo — 5soXo — llrfipoxo + Xq 
+ 13m2xQ + SsqXq + Srn^soXg) ^ 



(so - 2x0 + 1)2 

32>\]^ / 9 9 9 99 9\ 

— - — {—3mg — 5so + 2mgSQ + 3xo + Gm^xo + 3soXo — Xq — Sm^Xg) 

\/9 2 2 2 22\\ R'^iCn ('^)) C*"^ 

— 32A2(?TT,3 + So + 2m3So — xo + 2mgXo + 3soXo — 3xo — bm^Xn) 

^ So - 2xo + 1 

+ ^Ai(so-x2)|^(|cf(s)|2 + |Cio| 

'^^^/2 2 2 2\l ^-^'^ I 

+ ——{—2mg — 2so — 4m3So + xo + Sm^xo + soxo + m^soxo)- 



-128,. 2 ^ --2 ,^e(C|ff(s))C7effl 

+ + go - xo - m^xo - % '7 , 

3 So - 2xo + 1 J 
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T2{so,Xo) 



1 I / 112Ai 



40(^ + {\cf{sr + \Cio\') + (-64 + 160(^ 
(-1 + xo) + 16A2(-3 + 5xo)^ {\cf{s)\^ + \Cw 



+ A2) (1 + m 



c2\ 



\c: 



eff|2 



So - 2xo + 1 



+ 



/ 448Ai 



l-xo)+64A2(5xo-l) 



ir'eff|2 



So - 2xo + 1 



64X2Re{cfis))cf\ 



4) 



^ i\cf( 



+ \c- 



10 1 



+ 



-256 



(1 + m 



\C. 



eff|2 



So -2x0 + 1 



Here, x = sq — m^ + ie, Ai = Ai/m^ and A2 = ^2/^^. As the structure function T3 does not contribute 

rr 

to the branching ratio, we did not consider it in our present work. The Wilson coefficient Cg (s) 
depends both on the variables xo and sq arising from the matrix element of the Four-Fermi-operators. 

The branching ratio for B Xsi^i~ is usually expressed in terms of the measured semileptonic 
branching ratio Bgi for the decays B — > Xcii'i- This fixes the normalization constant Bq to be, 



3a' \V^sytb\' 



16vr2 \V,,f /( 



mclKimc) 



where 



f{mc 



1 



24 mt Inm-c 



(43) 



(44) 



is the phase space factor for T(B Xcii^e) and the function K{rhc) accounts for both the 0{as) QCD 
correction to the semileptonic decay width and the leading order (l/mfo)^ power correction [p^]. 
It reads as: 

K(mc) = IH g{mc) + ^ o , (45) 



2ml 



where 



9{rhc) 
Hrhc) 



f{rhc) 



f{rhc 



-9 + 24m^ - 72m^ + 72m^ - 15m^ - 72m^ Inm^ 



(46) 
(47) 



and the analytic form of j4o(m,c) can be seen in [Rl]. Note that the frequently used approximation 



|((vr'-f)(l 



' + 1) holds within 1.4% accuracy in the range 0.2 < z < 0.4. The equation 



g{z) = -1.671 + 2.04(z - 0.3) - 2.15(z - 0.3)^ is accurate for 0.2 < z < 0.4 to better than one per 
mille accuracy and that is what we have used here. 

The double differential ratio given in eq. (^) agrees in the (V — A) limit with the corresponding 
expression derived for the semileptonic decay B 
amounts to the following transcription: 



Xcii'e in |33] (their eq. (3.2)). Taking this limit 



a 



eff 



-Cin — — 
iu 2 



Gp a 
V2Tr 



eff 



Vcb 



(48) 
(49) 
(50) 



(42) 



14 



The hadron energy spectrum can now be obtained by integrating over sq. The imaginary part can 
be obtained using the relation: 



1 



(51) 



The kinematic boundaries are given as: 



max{'ml, —1 + 2xo + ^mf) < sq < Xq , 

1/ 2 

ms< xo < -(1 + 



4mf) 



(52) 



Here we keep rhi as a regulator wherever it is necessary and abbreviate Cg^ = Cg^(s = 1 — 2xo + rn^). 
Including the leading power corrections, the hadron energy spectrum in the decay B — > Xsi~^i~is 
given below: 



dxQ 



eff|2 , \r< 1 2 
9 I + l^lOl 



+ 



50 + -^Ifi'l + ^2^2 



eff|2 



xq - i(l +rTi2) 



,(7,9) 



.eff ^ ^eff 



dsn 



+ Aifcl^ 



(7,9). di^e(C9^ff) eff I ^ ,(7,9)d^^e(C|^ effl 
2 ' dso ' d4 I 



+ 6{xo - ^(1 + m2 _ 4m2))/5(Ai, A2) + <5'(xo - ^(1 + - 4m2))/5,(Ai, A2) 



(53) 



The functions g^^'^^^ , g^^ , Qi'^^ , /i^^^ , ^1'^'^^ , ^1^^ , ^j"^'^^ in the above expression are the coefficients of the 
power expansion for different combinations of Wilson coefficients, with g^^'^^ being the lowest 
order (parton model) functions. They are functions of the variables and rhs and are given in 
appendix |A|. The singular functions 5' have support only at the lowest order end point of the 
spectrum, i.e., at Xq"*^ = i(l + ml — 4ifif). The auxiliary functions fsiXi, A2) and /^'(Ai, A2) vanish 
in the limit Ai = A2 = 0. They are given in appendix ^ The derivatives of Cg^ are defined as 



eff 



eff 



-(s = 1 — 2xo + So; ^0 = "^s) {n = 1,2). In the {V — A) limit our eq. ( |53| ) for the hadron 



energy spectrum in ^ Xsi^l agrees with the corresponding spectrum in B ^ Xii^i given in [^] 
(their eq. (Al)). Integrating also over xq the resulting total width for B X^^^^^agrees again in the 



{V — A) limit with the well known result [16| 



The power-corrected hadron energy spectrum 



dB{B^Xse+f-) 
dEo 



(with Eq = nihXQ) is displayed in 



Fig. |5| through the solid curve, however, without the singular 6, 6' terms. Note that before reaching 
the kinematic lower end point, the power-corrected spectrum becomes negative, as a result of the A2 
term. This behavior is analogous to what has already been reported for the dilepton mass spectrum 



dB(B^Xse+e-) ■ 

d^ ' 



in the high region |jT^, signaling a breakdown of the ^ expansion in this region. The 
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Figure 5: Hadron energy spectrum '^^^^'d^p^ ^ parton model (dotted line) and including 

leading power corrections (solid line). Formi,/2 < Eq < rrih the distributions coincide. The parameters 
used for this plot are the central values given in Table^ and the default values of the HQET parameters 
specified in text. 

terms with the derivatives of Cg in eq. (^) give rise to a singularity in the hadron energy spectrum 
at the charm threshold due to the cusp in the function Y{s), when approached from either side. The 
hadron energy spectrum for the parton model is also shown in Fig. |5|, which is finite for all ranges of 

Eq. 

What is the region of validity of the hadron energy spectrum derived in HQET? It is known that 

in i? ^ Xsi^i^ decay there are resonances present, from which the known six Q populate the xq 

(or Eq) range between the lower end point and the charm threshold. Taking this into account and 

what has been remarked earlier, one concludes that the HQET spectrum cannot be used near the 

resonances, near the charm threshold and around the lower endpoint. Excluding these regions, the 

spectrum calculated in HQET is close to the (partonic) perturbative spectrum as the power corrections 

are shown to be small. The authors of ref. |^0[, [| who have performed an l/rUc expansion for the 

dilepton mass spectrum ^ ^ and who also found a charm-threshold singularity, expect a 

reliable prediction of the spectrum for < 3m^ corresponding to Eq > ^{1 -\-rh1 — 3m^) ~ 1.8 GeV. 

In this region, the effect of the l/rub power corrections on the energy spectrum is small and various 

spectra in B ^ X^^^^^ calculated here and in ref. |15] can be compared with data. 

^The 0{l/m^) correction to ^ has also been calculated in ref. |2l], however, the result differs in sign 

from the one in ref. |po| . It seems that this controversy has been settled in favor of ref. pc[ |. 
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The leading power corrections to the invariant mass spectrum is found by integrating eq. (41) with 



respect to xq. We have already discussed the non-trivial hadronic invariant mass spectrum which 
results from the 0{as) bremsstrahlung and its Sudakov-improved version. Since we have consistently 
dropped everywhere terms of 0{Xias) (see eq. (p|)), this is the only contribution to the invariant mass 
spectrum also in HQET away from sq = rh'^, as the result of integrating the terms involving power 
corrections in eq. (pl|) over xq is a singular function with support only at sq = m^. Of course, these 
corrections contribute to the normalization (i.e., branching ratio) but leave the perturbative spectrum 
intact for sq / m^. 

5 Hadronic Moments in 5 ^ X^^+f in HQET 

We start with the derivation of the lowest spectral moments in the decay B Xsi'^i~at the parton 
level. These moments are worked out by taking into account the two types of corrections discussed 
earlier, namely the leading power 1/mb and the perturbative 0{as) corrections. To that end, we 
define: 

f (^^0 - ^')"^^ T^f- d^od^o ' (^^) 
' ' Bq J dsodxo 

for integers n and m. These moments are related to the corresponding moments (xq^(so — rn^)"^) 
obtained at the parton level by a scaling factor which yields the corrected branching ratio B = 
BoM'f+f. Thus, 

{x^{so-^'sr) = ^M^-f. (55) 

The correction factor Bq/B is given a little later. We remind that one has to Taylor expand it in terms 
of the 0{as) and power corrections. The moments can be expressed as double expansion in 0{as) 
and 1/mi, and to the accuracy of our calculations can be represented in the following form: 

M\lf = D^^'"^^ + ^Cg^yl^"-'-) + XiDt'""^ + AsZ)^"''"^ , (56) 

with a further decomposition into pieces from different Wilson coefficients for i = 0, 1, 2: 

The terms j^"''^^ and S^^''^'^ in eq. ( [5^ ) result from the terms proportional to Re{Cg^)CY^ and |Cg^p 
in eq. (|4l|), respectively. The results for q,("'™)^^("''")^-)/|"'™')^(^("'™) are presented in appendix]^. Out 
of these, the functions a-"'"*^ and are given analytically, but the other two and S^^'"^^ are 

given in terms of a one-dimensional integral over xq, as these latter functions involve the coefficient 
which is a complicated function of xq. 
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The leading perturbative contributions for the hadronic invariant mass and hadron energy moments 
can be obtained analytically by integrating eq. (pO|) and eq. (p9[), respectively, yielding 



^(0,0) ^ ^5 - AtT^ ^(1^0) ^ ^(2,0) ^ _^ 

9 ' 675 ' 486 ' 

^(0,1) ^ 1381 - 210^2 ^(0,2) ^ 2257 - 320^^ 

1350 ' 5400 ■ ^ ^ 

The zeroth moment n = m = is needed for the normalization and we recall that the result for A^^'^^ 
was derived by Cabibbo and Maiani in the context of the 0{as) correction to the semileptonic decay 
rate B Xiu^ quite some time ago Likewise, the first mixed moment A^^'^^ can be extracted 
from the results given in |32] for the decay B Xlvi after changing the normalization, 



A(i'i) = -. (59) 
50 ^ ' 

For the lowest order parton model contribution D^'^\ we find, in agreement with |^2|, that the first 
two hadronic invariant mass moments {sQ — rfiD, {{sQ — m^)'^) and the first mixed moment {xQ{sQ — m1)) 
vanish: 

=0 forn = 1, 2 and D^^'^^ = . (60) 

We remark that we have included the s-quark mass dependence in the leading term and in the power 
corrections, but omitted it throughout our work in the calculation of the explicit ag term. All the 
expressions derived here for the moments agree in the V — A limit (and with rhg = in the perturbative 
as correction term) with the corresponding expressions given in |^2|. From here the full ©(a^rris) 
expressions can be inferred after adjusting the normalization Fq Bq'^C^. We have checked that a 
finite s-quark mass effects the values of the A^"'™-* given in eq. (58-59) by less than 8% for nis = 0.2 
GeV. 

We can eliminate the hidden dependence on the non-perturbative parameters resulting from the 
6-quark mass in the moments A^^+J™^ with the help of the HQET mass relation. As nis is of order 



^QCD, to be consistent we keep only terms up to order mj/mf^ |43|. An additional mfe-dependence 



is in the mass ratios 1111 = ^. Substituting nib by the B meson mass using the HQET relation 
introduces additional 0(l/mB,l/m^) terms in the Taylor expansion of eq. (|55[). We get for the 

^(0,0) 

^i+i- 



following normalization factor for B/Bq = Aif*'^^ " 



+ 



/ dxo—^{-^ml + %m\xl + Qmlxl - ^m\xl)\Cf^\^ 
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+ 



TT 3 niB 3 



+ 



B 



16 



C7- 



10 



dxo(64x^i?e(C|ff)C7«ff + ^(3 - 4xo)xl\cf\'') 



Ai 



+ 



9 J*-^7 + ~^ 

(4 + 91n(4^))C7«ff'-3C72o 
/' dxo(64(-l - 4x0 + 7xg)i?e(C|ff)Cfff + 16(-1 + 15xg - 20xg)|C|^|2 

JO 



(61) 



m 



B 



Here, the and terms result from the expansion of ln(4m^/m^). The first two moments and the 



rriB 



first mixed moment, {xq)B/Bq, {xq)B/Bq, {sq — 'mj)B/Bo, {{sq — rhf) )B /Bq and {xq{sq — 'mj))B/Bo 
are presented in appendix 

With this we obtain the moments for the physical quantities valid up to ©(a^/m^, 1/m^), where 
the second equation corresponds to a further use of = 0{Aqcd)- We get for the first two hadronic 
invariant mass moments ^ 



{Sh) = m2 + A2 + (m|-2AmB)(so-m2) + (2AmB-2A2-Ai-3A2)(xo), 

(Sff) = + 21^771^ + 2m^(m|-2AmB)(so-m^) + 2m^(2Amij-2A2-Ai-3A2)(xo) 

+ {m% - 4Am|)((so - rh^.f) + 4A''ml{xl) + 4Am|(xo(so - m^)) , (62) 
= (m|j - 4Am|)((so - m^?) + 4.A'^ml{xl) + 4Am|(xo(so - m^)) , 



and for the hadron energy moments: 



(Eh) 
{El) 



- A1+3A2 , ( , Ai + 3A2 ^ , , 

A^ + (2AmB - 2A2 - Ai - 3A2)(xo) 
+(m| - 2Ami3 + A^ + Ai + 3A2)(2;g) . 



(63) 



One sees that there are linear power corrections, ©(A/m^), present in all these hadronic quantities 
except {S'jj) which starts in 



cis A 
rriB ' 



5.1 Numerical Estimates of the Hadronic Moments in HQET 

Using the expressions for the HQET moments given in appendix |D|, we present the numerical results 
for the hadronic moments in i? ^ Xgl^C-" , valid up to 0{as/m\, 1/m^). We find: 



(xo) = 0.367(1 +0.148— -0.204 ^ 

vr TUB vr 



0.030- 



A 

niB 



A2 Ai A2 

0.017^ + 0.884^+3.652 



m 



B 



m 



B 



ni 



2 / ' 
B 



^Our first expression for (Sh), eq. (^), does not agree in tlie coefficient of (sq — ml) witli ttie one given in H] (tlieir 
eq. (4.1)). We point out that m% sfiould have been replaced by ml in this expression. This has been confirmed by Adam 
Falk (private communication). Dropping the higher order terms given in their expressions, the hadronic moments in 
HQET derived here and in |32] agree. 
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(4) 


a, A CKo A 
= 0.147(1 +0.324 ' 0.221 0.058 

vr mB TT rriB 


0.034^ + 1.206^ + 4.680^) , 


{xo{so - 




= 0.041 — (1 + 0.083 +0.124^+0.172^ 

vr niB m'^ 


! 


{so 


-ml) 


oi AAA 
= 0.093— (1 + 0.083 ) + 0.641^ +0.589^ 

TT niB TaB rriB 


! 






= 0.0071— (1 + 0.083— ) -0.196^. 

vr TUB 


(64) 



As already discussed earlier, the normalizing factor B/Bq is also expanded in a Taylor series. Thus, 
in deriving the above results, we have used 

B , Os A A^ Ai A2 , 

— = 25.277 (1 - 1.108— -0.083 0.041^ + 0.546^- 3.439^). 

Bq tt niB m^B "^B 

The parameters used in arriving at the numerical coefficients are given in Table | and Table |. 

Inserting the expressions for the moments calculated at the partonic level into eq. ( p2| ) and 
eq. (p3D, we find the following expressions for the short-distance hadronic moments, valid up to 
0(a,/m|,l/m|): 

{Sh) = m|(^ +0.093— -0.069— — +0.735— + 0.243^ + 0.273^ -0.513^), 

vr niB TT ms tti'^ 

{S%) = m|5(0.0071— + 0.138— — +0.587^ -0.196^), (65) 
vr niB TT 

(Eh) = 0.367mB(l + 0.148— -0.352— — + 1.691— + 0.012^5- + 0.024^ + 1.070^), 

vr niB TT ruB rn-^ rw^ 

,„n, , tts A a, A A^ Ai A2 , 

{E^h) = 0.147m|(l + 0.324— -0.128 ^ + 2.954 + 2.740^- 0.299^ + 0.162^). 

vr niB TT niB mB rn'^ 

Setting ms = changes the numerical value of the coefficients in the expansion given above (in which 
we already neglected a^mg) by at most 1%. With the help of the expressions given above, we have 
calculated numerically the hadronic moments in HQET for the decay B Xgi'^i^, £ = fj,,e and 
have estimated the errors by varying the parameters within their ±la ranges given in Table |. They 
are presented in Table | where we have used A = 0.39 GeV, Ai = -0.2 GeV^ and A2 = 0.12 GeV^ 
Further, using as{mf,) = 0.21, the explicit dependence of the hadronic moments given in eq. ( |65[ ) on 
the HQET parameters Ai and A can be worked out: 

A A^ A 

{Sh) = 0.0055m|(l + 132.61 +44.14^+49.66^), 

mB m^ 

A A^ A 

(Sl) = 0.00048m|}(l + 19.41 + 1223.41^ - 408.39^) , (66) 

mB 

A 

{Eh) = 0.3727715(1 + 1. 
{Efj) = 0.150m|(l + 2. 

Ill 
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A 


A2 


Ai 

+ 0.02^ 




+ 0.01^ 


mB 


mj^ 


m% 


A 


A2 


-0.29 A 




+ 2.68^ 


mB 




m% 



While interpreting these numbers, one should bear in mind that there are two comparable expansion 
parameters A/m^ and Os/vr and we have fixed the latter in showing the numbers. As expected, the 
dependence of the energy moments (-E^) on A and Ai is very weak. The correlations on the HQET 
parameters Ai and A which follow from (assumed) fixed values of the hadronic invariant mass moments 
(Sff) and (Sjj) are shown in Fig. We have taken the values for the decay B XsfJ.^fJ.~ from Table 
^ for the sake of illustration and have also shown the presently irreducible theoretical errors on these 
moments following from the input parameters mt, ctg and the scale /i, given in Table |l[ The errors 
were calculated by varying these parameters in the indicated range, one at a time, and adding the 
individual errors in quadrature. This exercise has to be repeated with real data in B ^ Xgi^i^to 
draw any quantitative conclusions. 

The theoretical stability of the moments has to be checked against higher order corrections and the 
error estimates presented here will have to be improved. The "BLM-enhanced" two-loop corrections 
p^ ] proportional to a^Po, where Pq = 11 — 2nj/3 is the first term in the QCD beta function, can be 
included at the parton level as has been done in other decays |32.45], but not being crucial to our point 
we have not done this. More importantly, higher order corrections in as and 1/m^ are not included 
here. While we do not think that the higher orders in will have a significant influence, the second 
moment {Sj^) is susceptible to the presence of l/m^ corrections as shown for the decay B Xiui 
[^ ]. This will considerably enlarge the theoretical error represented by the dashed band for {Sfj) in 
Fig. ^. Fortunately, the coefficient of the A/niB term in {Sh) is large. Hence, a good measurement of 
this moment alone constrains A effectively. Of course, the utility of the hadronic moments calculated 
above is only in conjunction with the experimental cuts. Since the optimal experimental cuts in 
B ATg^'''^" remain to be defined, we hope to return to this and related issue of doing an improved 
theoretical error estimate in a future publication. 

Related issues in other decays have been studied in literature. The classification of the operators 
contributing in 0{l/m^), estimates of their matrix elements, and effects on the decay rates and spectra 
in the decays B — > Xli^i and B (D, D*)ii'i have been studied in refs. |47-|49|]. Spectral moments 
of the photon energy in the decay B —>■ Xgj have been studied in ref. [^]. For studies of 0{l/ml) 
contributions in this decay and the effects of the experimental cut (on the photon energy) on the 



photon energy moments, see ref. [51|. 

Finally, concerning the power corrections related to the cc loop in B ^ Xgi'^l", it has been 
suggested in pO[| that an ©(Ag^^/m^) expansion in the context of HQET can be carried out to 
take into account such effects in the invariant mass spectrum away from the resonances. Using the 
expressions (obtained with mg = 0) for the l/m'^ amplitude, we have calculated the partonic energy 
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Figure 6: {Sh) (solid bands) and {Sjj) (dashed bands) correlation in (Xi-A) space for fixed values 
(Sh) = 1-64 GeV^ and (Sjj) = 4.48 GeV^ , corresponding to the central values in Table The curves 
are forced to meet at the point Ai = —0.2 GeV^ and A = 0.39 GeV. 



HQET 


(Sh) 




{Eh) 


{E%) 




(GeV^) 


(GeV^) 


(GeV) 


(GeV2) 


e'^e~ 


1.64 ± 0.06 
1.79 ±0.07 


4.48 ± 0.29 
4.98 ±0.29 


2.21 ±0.04 
2.41 ± 0.06 


5.14 ±0.16 
6.09 ±0.29 



Table 3: Hadronic spectral moments for B Xg^^ and B X^e^e^ in HQET with A = 0.39 GeV , 
Ai = —Q.2GeV^ , and A2 = 0.12GeV^ . The quoted errors result from varying fj.,as and the top mass 
within the ranges given in Table 



moments A(xq), which correct the short-distance result at order X^lfn 



B _ 256C72A2 /•V2(i-4mf) 
~Bq ~ 27m2 
1-2x0 
4m? ' 



dxoXn^^Re 



'Fir) {cfi?-2.,)^2cf--^±^^ 



(67) 



F(r) 



^Jr{\-r) 
1 



arctan 



1 — r 



in , ± ITT 



1 



< r < 1 , 
r > 1 . 



(68) 



[ 2v/r(r 1 + ^1- 1/r 

The invariant mass and mixed moments give zero contribution in the order we are working, with 
rUg = 0. Thus, the correction to the hadronic mass moments are vanishing, if we further neglect terms 
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proportional to and ^Aj, with i = 1,2. For the hadron energy moments we obtain numerically 



^{EH)i/ml = mBA(xo) = -0.007 GeV , 

)l/ml 



^{E%)i/m2 = m|A(x2) = -0.013 GeV2 , (69) 



leading to a correction of order —0.3% to the short-distance values presented in Table ^. The power 
corrections presented here in the hadron spectrum and hadronic spectral moments in B ^ Xgi'^i'sne 
the first results in this decay. 

6 Hadron Spectra and Moments in the Fermi Motion Model 

In this section, we study the non-perturbative effects associated with the bound state nature of the 
B hadron on the hadronic invariant mass and hadron energy distributions in the decay B Xsi'^i~ ■ 
These effects are studied in the FM model The hadronic moments in this model are compared 
with the ones calculated in the HQET approach for identical values of the equivalent parameters. We 
also define this equivalence and illustrate this numerically for some values of the FM model parameters 
resulting from fits of data in other B decays. With the help of the phenomenological profiles in the 



FM model, we study the effects of the experimental cuts used by the CLEO collaboration [30| on the 
hadron spectra and spectral moments in the decay B Xsi~^i~- The resulting branching ratios and 
the hadronic invariant mass moments are calculated for several values of the FM parameters and can 
be compared directly with data when it becomes available. 



6.1 Hadron spectra in i? — > Xsi~^i in the Fermi motion model pS 



The Fermi motion model [28| has received a lot of phenomenological attention in B decays, partly 



boosted by studies in the context of HQET showing that this model can be made to mimic the effects 
associated with the HQET parameters A and Ai p9| , p^ . We further quantify this correspondence in 
this paper. In the context of rare B decays, this model has been employed to calculate the energy 



spectra in the decay B Xg + 7 in [52|, which was used subsequently by the CLEO collaboration in 
their successful search of this decay [^]. It has also been used in calculating the dilepton invariant 
mass spectrum and FB asymmetry in — > Xsl^i^ in ref. |15]. 



The FM model has two parameters pp and the spectator quark mass niq. Energy- momentum 
conservation requires the 6-quark mass to be a momentum-dependent parameter determined by the 
constraint: 



ruhip) = mB +mq - 2mB\Jp'^ + ruq'^ ; p = \p\ ■ (70) 
The 6-quark momentum p is assumed to have a Gaussian distribution, denoted by </>(p), which is 

23 



determined hy pp 

4 -p2 

with the normahzation dpp^ (f){jp) = 1. In this model, the HQET parameters are calculable in 
terms oi pp and mq with 

A = J dpp'^(t){p)ym'^ + p"^, 

POO Q 

Ai = - / dppUip) = -ttPf • (72) 
Jo ^ 

In addition, for niq = 0, one can show that A = 2pf/\/tt. There is, however, no parameter in the FM 
model analogous to A2 in HQET. Curiously, much of the HQET malaise in describing the spectra in 
the end-point regions is related to A2, as also shown in For subsequent use in working out the 

normalization (decay widths) in the FM model, we also define an effective 6-quark mass by 



m. 



eff . 

b 



(/ dpp'mkipfmf''' ■ (73) 







The relation between m^, m;,. A, Ai and A2 in HQET has already been stated. With the quantity 
defined in eq. (73) and the relations in eqs. (72) for Ai and A, the relation 

mB = mf + K-\^/{2mf) , (74) 

is found to be satisfied in the FM model to a high accuracy (better than 0.7%), which is shown in 
Table ^ for some representative values of the HQET parameters and their FM model equivalents. We 
shall use the HQET parameters A and Ai to characterize also the FM model parameters, with the 
relations given in eqs. ( |72[ ) and ( [7^ ) and in Table ^. 

With this we turn to discuss the hadron energy spectrum in the decay B Xgi^i" in the FM 
model including the 0{as) QCD corrections. The spectrum -^^{B Xgi'^i^) is composed of a 
Sudakov improved piece from Cg and the remaining lowest order contribution. The latter is based on 
the parton model distribution, which is well known and given below for the sake of completeness: 



lOl^ 



T = So^\Umt-2m',ml + mt + mls + m',s-2s')(\cfis)\' + \C_ 

0.5 ^ 

16 IC^^P 
3 s 
+ 16{mt - 2mlml + - mis - mls)Re{cf{s)) C^^] , (75) 

J {ml + s — mIY — Am^s 
u = 2 ' 

BsiG^\V^Vul3a' , 
^0 = TTT^Zs — T^'^b' 



si 



1927r-^ 167r^ 



Tsi = ^l^^^h i^M^c). (76) 
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Note that in the lowest order expression just given, we have |Cg = |y(s)p + 2CQRe(Y{s)) with 

the rest of C^^{s) now included in the Sudakov- improved piece as can be seen in eq. (p2). To be 
consistent, the total semileptonic width Tgi, which enters via the normalization constant Bq, has also 
to be calculated in the FM model with the same set of the model parameters. We implement the 
correction in the decay width by replacing the 5-quark mass in given in eq. ( [76| ) by m^^. (See [^] 
for further quantitative discussions of this point on the branching ratio for the decay B — > Xgi^i^.) 
The hadronic invariant mass spectrum in the decay B — > Xsi~^i~iJi this model is calculated very much 
along the same lines. The kinematically allowed ranges for the distributions are mx < Eh < thb 
and mj^ < Sh < ^'b^ ^^'^ recall here that the physical threshold has been implemented by 
demanding that the lowest hadronic invariant mass produced in the decay B X<j^"'"£~ satisfies 
mx = max{mK,mq + 'ms)- The results for the hadron energy and the hadronic invariant mass spectra 
are presented in Figs. ^ and |9|, respectively. We do not show the Sh distribution in the entire range, 
as it tends monotonically to zero for larger values of Sh- 

' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' 



if 




Eh [GeV] 

Figure 7: Hadron energy spectrum in B ^ Xsi'^i~in the Fermi motion model based on the per- 
turbative contribution only. The solid, dotted, dashed curve corresponds to the parameters (Ai,A) = 
(-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV^, GeV), respectively. 

A number of remarks is in order: 

• The hadron energy spectrum m. B ^ Xsl'^l~''\.s rather insensitive to the model parameters. Also, 
the difference between the spectra in the FM and the parton model is rather small as can be 
seen in Fig. 0. Since, away from the lower end-point and the cc threshold, the parton model 
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m 

CD 




Figure 8: Hadron energy spectrum in B ^ Xgl I based on the perturbative contribution only, in 

eff 

the Fermi motion model (dotted curve) for {pp,mq) = (252,300) {MeV, MeV), yielding = 4.85 
GeV, and in the parton model (long-short dashed curve) for ruh = 4.85 GeV. 

and HQET have very similar spectra (see Fig. the estimates presented in Fig. ^ provide a 
good phenomenological profile of this spectrum for the short-distance contribution. Very similar 
conclusions were drawn in |^3| for the corresponding spectrum in the decay B Xuii^e, where, 
of course, the added complication of the cc threshold is not present. 

• In contrast to the hadron energy spectrum, the hadronic invariant mass spectrum in B ^ 
Xsl^£~ is sensitive to the model parameters, as can be seen in Fig. |^. Again, one sees a close 
parallel in the hadronic invariant mass spectra in B ^ X^^+^^and B X^ii^i, with the latter 
worked out in |34|. We think that the present theoretical dispersion on the hadron spectra in 
the decay B X^^^^^can be considerably reduced by the analysis of data in B ^ X^ii^e- 

• The hadronic invariant mass distribution obtained by the 0(as )-corrected partonic spectrum 
and the HQET mass relation can only be calculated over a limited range of Sh, Sh > itt-b^, as 
shown in Fig. ^. The larger is the value of A, the smaller is this region. Also, in the range where 
it can be calculated, it depends on the non-perturbative parameter nih (or A). A comparison of 

rr 

this distribution and the one in the FM model may be made for the same values of rub and . 
This is shown for mi, = 4.85 GeV in Fig. ^for HQET (long-short dashed curve) to be compared 
with the dotted curve in the FM model, which corresponds to = 4.85 GeV. We see that 
the two distributions differ though they are qualitatively similar. 
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Figure 9: Hadronic invariant mass spectrum in the Fermi motion model and parton model, based 
on the perturhative contribution only. The solid, dotted, dashed curve corresponds to the parameters 
(Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV^ , GeV), respectively. The parton model (long- 
short dashed) curve is drawn for rub = 4.85 GeV. 

6.2 Numerical Estimates of the Hadronic Moments in FM model and HQET 

To underline the similarity of the HQET and FM descriptions in B ^ Xgi'^i^, and also to make 
comparison with data when it becomes available with the FM model, we have calculated the hadronic 
moments in the FM model using the spectra just described. The moments are defined as usual: 

(X^) = ( / X-^^dXH)/B for X = S,E. (77) 

The values of the moments in both the HQET approach and the FM for n = 1, 2 are shown in Table |5| 
for the decay B Xgfi^fi" , with the numbers in the parentheses corresponding to the former. They 
are based on using the central values of the parameters given in Table || and are calculated for the 
same values of the HQET parameters A and Ai, using the transcriptions given in eqs. ([7^). Both the 
HQET and the FM model lead to strikingly similar results for the hadronic moments shown in this 
table. With {Sh) — (1.5 — 2.1) GeV, the hadronic invariant mass spectra in i? — > Xsi'^i~are expected 
to be dominated by multi-body states. 



27 



PF^rUq (MeVjMeV j 




Ai (GeV ) 


A (GeVj 


(450,0) 


4.76 


-0.304 


0.507 


(252,300) 


4.85 


-0.095 


0.422 


(310,0) 


4.92 


-0.144 


0.350 


(450, 150) 


4.73 


-0.304 


0.534 


(500, 150) 


4.68 


-0.375 


0.588 


(570, 150) 


4.60 


-0.487 


0.664 



eff 

Table 4: Values of non perturbative parameters , Ai and A for different sets of the FM model 
parameters {pF,mg) taken from various fits of the data on B Xs + (J/ip,^) decays discussed in 
ref 





{Sh) 




{Eh) 


{ED 


(Ai,A) in (GeV2, GeV) 


(GeV^) 


(GeV^) 


(GeV) 


(GeV2) 


(-0.3,0.5) 
(-0.1,0.4) 
(-0.14,0.35) 


2.03 (2.09) 
1.75 (1.80) 
1.54 (1.49) 


6.43 (6.93) 
4.04 (4.38) 
3.65 (3.64) 


2.23 (2.28) 
2.21 (2.22) 
2.15 (2.18) 


5.27 (5.46) 
5.19 (5.23) 
4.94 (5.04) 



Table 5: Hadronic spectral moments for B XsP-^n in the Fermi motion model (HQET) for the 
indicated values of the parameters (Ai,A). 

7 Branching Ratios and Hadron Spectra in B ^ Xgt^l^ with Cuts 
on Invariant Masses 

The short-distance (SD) contribution (electroweak penguins and boxes) is expected to be visible 
away from the resonance regions dominated by i? — > Xs{J/ip,ip' , ...) — > Xgl^i^. So, cuts on the 
invariant dilepton mass are imposed to get quantitative control over the long-distance (LD) resonant 



contribution. For example, the cuts imposed in the recent CLEO analysis |3C] given below are typical: 



cut A : < {nij/^ - 0.1 GeV)^ = 8.98 GeV^ , 
cut B : g2 < (mj/^ - 0.3 GeV)^ = 7.82 GeV^ , 
cut C : g2 > (m^/ + 0.1 GeYf = 14.33 GeV^ . 



(78) 



The cuts A and B have been chosen to take into account the QED radiative corrections as these 
effects are different in the e~^e~ and ^'^ ^~ modes. In a forthcoming paper [^], we shall compare the 
hadron spectra with and without the B {J/tp, •••) ~^ Xgi'^l" resonant parts after imposing these 
experimental cuts to quantify the theoretical uncertainty due to the residual LD-effects. Based on 
this study, we argue that the above cuts in greatly reduce the resonant part. Hence, the resulting 
distributions and moments with the above cuts essentially test (up to the non-perturbative aspects) 
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the SD contribution in B ^ Xs£^i~ . 

As mentioned in [^], the dominant BB background to the decay B — > ^^^"^^"comes from two 
semileptonic decays oi B oi D mesons, which produce the lepton pair with two undetected neutrinos. 
To suppress this BB background, it is required that the invariant mass of the final hadronic state 
is less than t = 1.8 GeV, which approximately equals mo- We define the survival probability of the 
B signal after the hadronic invariant mass cut: 

S{t) ^{1 -^dSH)/B , (79) 

and present S{t = 1.8 GeV) as the fraction of the branching ratio for B Xg^^^^ surviving these 
cuts in Table ^. We note that the effect of this cut alone is that between 83% to 92% of the signal for 
B Xsfi^n" and between 79% to 91% of the signal in B ^ Xse^e~ survives, depending on the FM 
model parameters. This shows that while this cut removes a good fraction of the BB background, 
it allows a very large fraction of the B X^^"''^" signal to survive. However, this cut does not 
discriminate between the SD- and LD- contributions, for which the cuts A - C are effective. 

With the cut A (B) imposed on the dimuon (dielectron) invariant mass, we find that between 57% 
to 65% (57% to 68%) of the B X^^^^^ signal survives the additional cut on the hadronic invariant 
mass for the SD contribution. The theoretical branching ratios for both the dielectron and dimuon 
cases, calculated using the central values in Table |l| are also given in Table ^ As estimated in |]l^] , 
the uncertainty on the branching ratios resulting from the errors on the parameters in Table || is 
about ±28% (for the dielectron mode) and ±21% (for the dimuon case). The wave- function-related 
uncertainty in the branching ratios is negligible, as can be seen in Table ^. This reflects that, like 
in HQET, the corrections to the decay rates for B — > Xsl'^(-~ and B Xlv^ are of order l/rn^, 
and a good part of these corrections cancel in the branching ratio for B Xgi^i". With the 
help of the theoretical branching ratios and the survival probability S{t = 1.8 GeV), calculated for 
three sets of the FM parameters, the branching ratios can be calculated for all six cases with the 
indicated cuts in Table |^ This gives a fair estimate of the theoretical uncertainties on the partially 
integrated branching ratios from the S-meson wave function effects. This table shows that with 10^ 
BB events, 0(70) dimuon and O(IOO) dielectron signal events from B X^^^^^should survive the 
CLEO cuts A (B) with m{Xs) < 1.8 GeV. With cut C, one expects an order of magnitude less 
events, making this region interesting for the LHC experiments. We show in Fig. |l^ hadron spectra 
in i3 — > Xs£~^i~, £^ = e^,iJ,^, resulting after imposing the CLEO cuts A, B,C, defined in eq. (IT^) . 
One sees that the general features of the (uncut) theoretical distributions remain largely intact: the 
hadron energy spectra are relatively insensitive to the FM parameters and the hadronic invariant 
mass spectra showing a sensitive dependence on them. Given enough data, one can compare the 



experimental distributions in B Xgt £ directly with the ones presented in Fig. IC . 
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FM parameters 
(Ai, A) in (GeV^, GeV) 


B ■ 10-6 


B ■ 10-6 


No s-cut 


No s-cut 


cut A 


cut B 


cut C 


cut C 

e"^e 


(-0.3,0.5) 


5.8 


8.6 


83% 


79 % 


57% 


57% 


6.4% 


4.5% 


(-0.1,0.4) 


5.7 


8.4 


93% 


91 % 


63% 


68% 


8.3% 


5.8% 


(-0.14,0.35) 


5.6 


8.3 


92% 


90 % 


65% 


67% 


7.9% 


5.5% 



Table 6: Branching ratios for B X^i'^l^ , i = ij,,e for different FM model parameters are given in 
the second and third columns. The values given in percentage in the fourth to ninth columns represent 
the survival probability S{t = 1.8 GeV), defined in eq. (fz^j, with no cut on the dilepton invariant mass 
and with cuts on this variable as defined in eq. 



We have calculated the first two moments of the hadronic invariant mass in the FM model by 
imposing a cut Sh < t^ with t = 1.8 GeV and an optional cut on q^ . 



{Sli) = ( 



d'^Bcutx ^ 2^ 



dSndq^ 



dSHdq')/{ 



dSudq^ 



dSudq^ 



for n = 1, 2 



(80) 



Here the subscript cutX indicates whether we evaluated {Sh) and {S'jj) with the cuts on the invariant 



dilepton mass as defined in eq. (78), or without any cut on the dilepton mass. The results are collected 
in Table |^ The moments given in Table | can be compared directly with the data to extract the 
FM model parameters. The entries in this table give a fairly good idea of what the effects of the 
experimental cuts on the corresponding moments in HQET will be, as the FM and HQET yield very 
similar moments for equivalent values of the parameters. The functional dependence of the hadronic 
moments on the HQET parameters taking into account the experimental cuts still remains to be 
worked out. 



FM 


No s-cut 


No s-cut 


cut A 


cut B 


cut C 


parameters 






e+ 


e 






e+ 


e 




(Ai,A) 


{Sh) 


{Si) 


{Sh) 


{SD 


{Sh) 


{Si) 


{Sh) 


{SD 


{Sh) 


{SD 


GeV2, GeV 


GeV2 


GeV^ 


GeV2 


GeV^ 


GeV2 


GeV^ 


GeV2 


GeV^ 


GeV2 


GeV^ 


(-0.3,0.5) 


1.47 


2.87 


1.52 


3.05 


1.62 


3.37 


1.66 


3.48 


0.74 


0.69 


(-0.1,0.4) 


1.57 


2.98 


1.69 


3.37 


1.80 


3.71 


1.88 


3.99 


0.74 


0.63 


(-0.14,0.35) 


1.31 


2.34 


1.38 


2.55 


1.47 


2.83 


1.52 


2.97 


0.66 


0.54 



Table 7: {Sh) o-nd {SD for B — > Xgi'^i , £ = fj,,e for different FM model parameters and a hadronic 
invariant mass cutSn < 3.24 GeV^ are given in the second to fifth columns. The values in the sixth to 
eleventh columns have additional cuts on the dilepton invariant mass spectrum as defined in eq. (f?^. 
The Sn-moments with cuts are defined in eq. ^^). 
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8 Summary and Concluding Remarks 



We summarize our results: 

• We have calculated the 0{as) perturbative QCD and leading 0{l/mi,) corrections to the hadron 
spectra in the decay B Xgi'^i' , including the Sudakov-improvements in the perturbative part. 

• We find that the hadronic invariant mass spectrum is calculable in HQET over a limited range 
Sh > TUB A and it depends sensitively on the parameter A (equivalently mh). These features 
are qualitatively very similar to the ones found for the hadronic invariant mass spectrum in the 



decay B XJi^e H] 



The 1/rn-fe-corrections to the parton model hadron energy spectrum in S — > Xsi~^i~ are small 
over most part of this spectrum. However, heavy quark expansion breaks down near the lower 
end-point of this spectrum and close to the cc threshold. The behavior in the former case has a 
similar origin as the breakdown of HQET near the high end-point in the dilepton invariant mass 
spectrum, found in ref. |15]. 



We have calculated the hadronic spectral moments (5*^) and {E^) for n = 1, 2 using HQET. The 
dependence of these moments on the HQET parameters is worked out numerically. In particular, 
the moments (5"^) are sensitive to the parameters A and Ai and they provide complementary 
constraints on them than the ones following from the analysis of the decay B — > Xivf^. The 
simultaneous fit of the data in i3 — > X^^+^^and B Xiv^ could then be used to determine 
these parameters very precisely. This has been illustrated in ref. pH] based on the present work. 



The corrections to the hadron energy moments l\{EH)i/ml ^'^'^ ^{^H)i/m'^ from the leading 
©(Ag^^^i/m^) power corrections have been worked out, using the results of ref. [^0|. We find 
that these corrections are very small. The corresponding corrections in A{S'^)i/^2 vanish in the 
theoretical accuracy we are working. 

We think that the quantitative knowledge of A and Ai from the moments can be used to remove 
much of the theoretical uncertainties in the partially integrated decay rates in i? — > Xuii^i and 
B Xsi'^i~ ■ Realating the two decay rates would enable a precise determination of the CKM 
matrix element Vub- 

As a phenomenological alternative to HQET, we have worked out the hadron spectra and spectral 
moments in B ^ Xgl'^l" in the Fermi motion model [^]. This complements the description of 
the final states \n B ^ presented in |T^, where the dilepton invariant mass spectrum 

and FB asymmetry were worked out in both the HQET and FM model approaches. We find 
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that the hadron energy spectrum is stable against the variation of the FM model parameters. 
However, the hadronic invariant mass is sensitive to the input parameters. For equivalent values 
of the FM and HQET parameters, the spectral moments are found to be remarkably close to 
each other. 

• We have worked out the hadron spectra and spectral moments in the FM model by imposing 
the CLEO experimental cuts designed to suppress the resonant cc contributions, as well as the 
dominant BB background leading to the final state BB —>■ Xsi~^£~ (+ missing energy). The 
parametric dependence of the resulting spectra is studied. In particular, the survival probability 
of the B Xst^i~ signal is estimated by imposing a cut on the hadronic invariant mass 
Sh < 3.24 GeV^ and on the dilepton invariant mass as used in the CLEO analysis. The spectra 
and moments can be directly compared with data. 

We hope that the work presented here will contribute to precise determinations of the HQET 
parameters and Vuh using the inclusive decays B Xst^i~ and B Xu£u£ in forthcoming B 
facilities. 
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Appendices 



A Coefficient Functions gf\ hf\ /.f'^^ /cf \ /cf''^ 

These functions enter in the derivation of the leading (l/m^) corrections to the hadron energy spectrum 
in i? ^ Xs('^(~ , given in eq. (|5^). 



gj^'"'^ = ^xl - mf^{-2rhl + 3xo + Sm^xo - Axl) , (A-1) 

(9,10) _ 1 fn^'^ I OQ^4 10^2^^ 10^2 c:o»5,2^2 , qr^S , on^^^N 



^va,iu; ^ / + ^^^^ ~ ISmfxQ - 18xg - 52m>g + 36xg + 20xg) , (A-2) 

Vq — ml 

^(9,10) ^ J^(3^2 _^ 23m^ _ - 21m^xo - 6xg - b2rhlxl + 36xg + 20x^) , (A-3) 

Xq — rhi 

= \J Xq — m^ — (lOm^ + lOrn^ — 3xo — ISm^xo — 3m^xo + 2xq + 2mgXQ) , (A-4) 
= ^ YT 7377^^(9^? + 34m^ + 104m^ + 110mf + 31mJ°-132mfxo-312m^xo 



- ISOmfxo - 18xg - 170m^xg - 58m^x^ + 74mf xg - 20mf xg + 72x1] + 564m^xg + 576m^xj] 

+ 22%rhlxl - llQxl - QlQmlxl - 436m^x^ - 20m^x^ + 72x^ + 240m^x^ + 24m^x^) , (A-5) 

1 1 16 
^ = , . (27m2 + 93m^ + 97m^ + 31m^ - 3xo - 63m^xo - 189m^xo 



129mf xo - 18xg - 108m^xg - 62mfx^ - 20m^x^ + 72xg + 324m^x[] + 180m^x[] - 60x^ 
152m^x^ - 20m^x^) , (A-6) 



= \lxl- rhll2S{-2rhl + xo + m^xo) , (A-7) 

^(7,9) ^ ^ ^ U{ml + 3mf + 2m^xo - 2x1 - ^^^s^g) , (A-8) 
Xq — ml 

g^''^^ = I =64(5m^ + 9mf — xq + 5rn^xo — 6xq — 12mgXQ) , (A-9) 
xl-ml 

= — 7xg-m2(-12m2-6m^ + 9xo + 19m^xo + 3xg + 15m^xg-28x|]) , (A-10) 
9 

= — ^xg-m2(-6mf + 3xo + 5m^xo + 3xg + 15m^xg-20xg) , (A-11) 



64 



Xq - 




Xq - 


m2(- 


y^o ■ 


— m^i 


y^o ■ 


— m^i 




3 



128 

h\ = -\/xq — m^(— 8m^ — 2m^ + 3xo — 3m^xo + 5xq + 5m^XQ) , (A-12) 

3 

h2'^^ = — r— \/xq — 7fi'j{—Am1 — 6m^ + 3xo — 15m^xo + 7xq + ISm^Xg) , (A-13) 



(9,10) 
1 

^(7,9) ^ a/x§ - m2 (-2mf + Xo + m-s^o) • (A-15) 



^(y.iuj ^ — Jxg-m2 (2m^-3xo-3m>o + 4x5) , (A-14) 
9 * 
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B Auxiliary Functions /^(Ai, A2), /,5'(Ai, A2) 



The auxiliary functions given below are the coefficients of the singular terms in the derivation of the 
leading (1/m^) corrections to the hadron energy spectrum in B ^ Xsi~^£~ , given in eq. 



-(l-rfi2)3(5-m2)Ai 



+ -(l-m2)3(_l + 5^2)^^ 



(|c. 



+ 



(1 + 12rhf - 88mf - Amj. - 36mfmi - 736rhfmi + 5m^ + 24mf + 720m^m 



eff|2 I |2 
9 I + |t-^io| 



,2-2 



4-2 



,2-4 



,4-4 



+ 24rnfrh^ + IGOrfifm^ — 5m^ — SQrhfrh^ — 56mfrfif)^r^ + -(—1 + 3 

TTlj O 



+ Umf - 2ml + 166mf + 8mf + 154mf + 2mf + 50mf - 5mf )A2 
+ 



I a 



eff|2 



mr 



-(1 - m2)3(7 + ^2)^^ ^ _ ^2)3(13 + 15m2)A2 



+ Ai(-l + mf)^(- 



.2d\cf\' , 8 di?e(C|ff) ^eff 



9 







~'~ 3 dsQ 



(B-16) 



/<s'(Ai, A2) 



1 



fioAi^-(l-m2)5(|c; 



efr|2 

9 I 



+ \Cio\^ 



2 \C' 



«ff|2 



mf 



(B-17) 



C The Functions a^, 7i, 



The functions entering in the definition of the hadron moments in eq. ( p7\ ) are given in this appendix. 
Note that the functions a^"'™"^ and Z?!"'™"^ multiply the Wilson coefficients jC^^p and C^q, respec- 
tively. Their results are given in a closed form. The functions ji"''"^^ multiply the Wilson coefficients 
Cf^Re{C^^), of which Re{C^^) is an implicit function of xq. Likewise, the functions si"^'"^^ mul- 
tiply the Wilson coefficient |Cg^p. The expressions for -)/|"''"^ and si"''^^ are given in the form of 
one-dimensional integrals over xq. 
The functions a,^"'"*^ 



a, 



(0,0) 



16 



32 



--(-8 - 26mf + 18mt + 22m° - llm°) + —{-1 + )^(1 + mf) ln(4m;^ 
y o 

64 

+ — m^(-9 - 2m^ + m^) In(ms) , 



(0,0) 



(0,0) 
' 



a, 



(0,0) 



(-4 + 38m2 - 42m^ - 26m^ - 15m^) + 16(-1 + mf)^(3 + 8mf 5rfi^) ln(4m;' 



^6 



^,2\2/ 



(C-18) 
(C-19) 
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+ 32rh1{-8 - 17ml " 2m^ + Hitt-s) , (C-20) 

4°'^) = -(-41-49m^ + 256mf-128m^-43m^) + — (-l + m^)3(l + m^)^ln(4rfi^) 
y o 

16 

+ -m;(3-m^2™j)l„(,n.). (C-21) 

(0,1) (0,0) 

a^°'^) = 1(21 + 167ml + 128m^ - 276m^ - 319mf) + y (-1 + m^)^(3 + 14m^ + 21m^ + lOrfi^) In 
32 

+ -m?(3-24™j-18,hj+™j)ln(m.), 
ai°'^^ = — (-119-144m^ + 45m^ + 320m^+45mf) + -(-l + m^)^ln(4m^) 

41:0 O 

- ymf(8 + 3m2)ln(m,) , 



aS°'^^ = 7^(-127 - 278m2 + 1075m^ - 800m^ + 49m^) + ^(1 - mlf{-7 - 17ml - + ^^t) H'^rnj 



27' * * * 9 

+ -mj{18 - 38ml - l^mf ) Hrris) , 
9 

4°'^) = 1(27 - i6ml + 1681m^ - 688m^ - 1189rnf ) + !(-! + mjfiS + 20m^ + 25m^) ln(4?7i^ 
9 3 

- ^m^(18 + 5Aml + 47mf ) ln(m5) , 
= 0, 

= -(-23 - 159m^-112m^ + 304m^-45m^)-— (-l + m^)^(l + m^)ln(4m^) 
y o 

+ — m5(— 39 — 7771^ + Grfig) In(ms) , 



^(1,0) ^ -(-93-469m2 + 704mf-127mf) + ^(-l + rn^)^(3 + 8?fi^ + 5m^)ln(4mf^ 
9 3 
112 

— m^(3 + 3rn^ + 2mj) In(ms) , 



= 0, 

aj^'^) = ^(-4 - 131m^ + 307TOf - 416m^ + 178m^) - -(-1 + ml)'^{l + 6ml + ^rhj) ln(4mf ) 
27 9 

+ -^771^(9 — 35ml ~ '''"^s) l^(™s) ) 

2 8 
4^'^) = -(-60 - 185m^ + 173m^ + 160m^ + 70mf ) + -(-1 + mlf{l + rhlf{3 + 5?fi^) ln(4rfif) 
y o 

+ -^m^(3 — 21m^ — 13m^) In(ms) , 

o oo 

13f 
64 



= —(119 - 176m^ - 1085m^ + 400m^ + 835mf ) + -^-(l - rhlf{l + ml) ln(4?fif) 
-Loo y 



— m°(28 + 5mj)ln(m,) , 



a(^'°) = 0. 



The functions pi'^'"^^ 



C-22) 



4mf) 



C-23) 



C-24) 



,2\ 



C-25) 

C-26) 
C-27) 

C-28) 

C-29) 
C-30) 

C-31) 

C-32) 
C-33) 

C-34) 
C-35) 



/3o°'°^ = ^(l-8m^ + 8m^-m«-24m^ln(m,)), (C-36) 
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_ 1 oim 

= -3 + 8m^ - 24mi + 24m^ - 5m^ - 24mf In(rfis) , 

= ^(7 - 25m^ + 160mf - 160m^ + 25m^ - 7mJ° + 120mj In(m^) + 120m^ In(m^)) , 

= 

= ^m^(7 - 20rn^ + 20m^ - 7mf + 24 In(rfis) - 48rfi^ In(ms)) , 
2 

= ^(2 - - 30mf + 30mf + 3mf - 2ml'^ - 120m^ 111(771^)) , 

= ;^(43 - 135m^ + 1260m^ - 1440m^ + 405mf - 153mJ° + 20™^^ + iQSOmj In(ms) 

+ 840rfi^ In(rfis)) , 

= ^(13 - 315m^ + 1500m^ - 1560rfi^ + 315mf + 147m^° - lOOm^^ + 350^4 i^^^^^ 



+ 840m° 111(7713)) , 

= 0, 

= ^(13 - 135m^ - 160mf + 320mf - 45m^ + 7mJ° - 600?fi^ In(ms) - 120mf In(ms)) , 

= ^(3 - 9771^ + 16771^ - 48777^ + 45771^ - 7777^° + 24777^ ln(7775) - 72777^ ln(77l3)) , 

= 0, 

= (23 - 45771^ + 1080771^ - 1440771^ + 5857nf - 24377i^° + 40774^^ + 10807n^ \n{ms) 

+ 600777^ ln(777s)) , 

= ^(13 + 45777^ - 120771^ - 45777^ + 147777^° - 40777^^ + 360777^ ln(77i5) - 600771^ \a{ms)) 



= 0, 



1 f\ 

— (-1 + 97n^ - 457n^ + 457nf - 97nf + ttiJ^ - 120771^ \n{ms)) , 

LOO 

. 



The functions ^y^"'"*^ 



Note that in the expressions given below Cg^ = Cq^(s = 1 — 2xo + 777^). The lower and upper 



limits of the xo-integrals are: xl^''"' = rhs and x™"^ = ^(1 + 777^ — 4r7i^) 



(0,0) 

7o 



max 

= 128/° dx^Jxl - ml{-2rhl + xq + mlxo)Re{C^^) , (C-54) 

7r = yr, (c-55) 



I ^ U4t //^9 Aid. a9 9 a99 

72 ' ' = / dxp — (4mg + lArrig — Xq — rrigXQ — 1277X5X0 — 4:Xq — 22771^X0 + 7xq 

I ^min / 2 ^9 



2 

(0,0) _ r°""^. ^ 

+ 157n^xg)ii:e(C|^) , (C-56) 

= 128/° dxQXoJxl - ml{-2ml + xq + mlxQ)Re{cf^) , (C-57) 



(0,1) 

7o 
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= ^r, (c-58) 

72°'^^ = — / dxQ — (4rng + Gm^ + Qm^xp + STmgXp — Sxp — lAm^XQ — STyfigXp 



y^o ~ 

- 9x|] - Slrhlxl + 28x^ + 60m^x^)i?e(C|^) , (C-59) 

max 

7^°'^^ = 128/° dxoxlJxl - m2(-2m^ + xq + rfi^xo)i?e(C|^) , (C-60) 

^(0,2) — dxoy Xq — m'^{—Amg — lOm^xo + 2mgXo + 6xq + IGinlxQ — 5x^ 

3 ixg^'" 

- 5m2x3)i?e(C|ff) , (C-61) 
72°'^^ = — / dxp — , ° (Srrtg + 12mg + Sm^xp + 72mgXo — 3xq — 25'mlxQ — TSm^Xg 

- Qxl - mmlxl + 35x^ + 75m^x^)i?e(C|^) , (C-62) 
7o^'°^ = 0, (C-63) 
7^'°^ = 128/" dxQ 7° ~ ' =(-2m^ + ?ngxp + mfxp + 2m^g;g - - mga;g)i?e(C|^) , (C-64) 

7^^'°) = ° dxpVx§-m2(-4m^-6m^ + 3xp-15m^xo + 7xg + 15m^x§)it!e(C'|^) , (C-65) 

3 Ja;^»" * 

7?''^ = 0, (C-66) 
128 r^o^^^ ]^ 

7^^'^^ = — — / dxo — (47?7,g + 4mgXp — 2r?7.gXp — 3rngXp — ITm^Xp + mgXp + 7mgXp 

•^^^ ' 'xi — mi 



+ 3x^ + ISm^x^ - - 5mlxl)Re{C^^) , (C-67) 

7^^'^^ = — - / ° dxQXoJxl - ml{-4.rhl - 6mf + 3xp - ISm^xp + 7x§ + lhrh1xl)Re{cf-) , (C-68) 
3 Jxg'*" 

7^^ = 0, (C-69) 

7^'°) = ° dxoJxl - m\{-2m\ + m^xp + m^xp + 2mlxl - xg - mlxl)Re{cf-) , (C-70) 
3 Jx^^"^ 

7^''°^ = 0. (C-71) 
The functions (^f '"'^ 

^(0,0) ^ 32 r^o ^^^f^Z~2(^_2rkl + Sxo + 3mlxo-4xl)\cf\\ (C-72) 

= ^^S"■^ (c-73) 

^(0,0) _ / — (6mg — xp — Qm^xp — 12mgXp + Bm^Xp + 15xp + ISm^Xp 

- 2Qxt)\cf\\ (C-74) 

5^'^'^ = /° dxQXQxIxl - mj{-2m1 + 3xp + 3rn^xo - 4xo)|C|^|^ , (C-75) 

3 Ja;^»" * 

5?^'^ = (C-76) 
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62''^^ = — / dxo — i^^^^={6m^ — Sm^xo + ISrhgXo — 3xq — SOrh^XQ — 57mgXQ + 3xq 



^(0,2) ^ 32 p /^2 _ ^2(_2^2 + 3^^ ^ 3^2^^ _ 4^2)|c.eff|2 ^ 

16 /"^cr"^ 



+ 20rt3)|C|Sr^ 



"^^0 ^xf^-mj 



dxQ ^L (-2m^ + 3m^a;o + 3m^a;o - 2m^x^ - 3a;^ - Sm^x^ + 4x^)|C|^|^ , 

1 — 



.(1,0) 


= 0, 




^(1,0) 


32 

y 


^max 


.(1,0) 
"2 


32 

y 


^max 


^(1,1) 


= 0, 






32 


^max 


y 


1 ^min 



dxQ^xl - rhl{-Qrn^g + 3xo + hrhlxQ + ^xl + Ihrhlxl - 20xo)|C|^|^ , 



I dxQ — . (4rng — Qm^gXQ — ^rh^gX^ — Ibm^x^ + 27inlxQ + 21ihgXQ 



+ 9x^ - dm'ix^ - 27x1 - ISm^a;^ + 20x^)|C|**|2 , 

(^2^'^^ = -;r dxQXoJ Xq — m'^{—6rhg + 3xo + brh^XQ + 3xo + ISm^XQ — 20xo)|C| 

3 Jii"™ 



= 0, 



^(2,0) _ f dxQxj Xq — ml{—2'm^ + 3?7i^xo + 3m^xo — 2fh?gX^ — 3^0 

9 7a;?r'*" 



3m?xg + 4x^)|cf|2 



4^'°) = 0. 



D Lowest Hadronic Moments (Parton Level) 

^^•^^^ = ^(-41m|-49m2-24(m|-m^)ln(4^))CfffV^(7m|-25mf)Ci^^ 



+ 



f\{\+ml/m%) 64 , n n n n 9 9. . pff. pff 

' dxQ^XQ{-mg - ArrigXQ + 2mBXQ + 2mgXQ)Re{Cf^)C^^ 

ms/mB f^B 

■■\{l-\-m1/m?g) 26 fF 

dxQ- — 2"a^o(— 3m^ + Qm\x^ + 6777,5X3 — 8m^Xo)|Cf ^ 



TT 3 3 miy 



^(l+31n(44))C?ff%^ 



mi 
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+ 

+ 
+ 



Jj dxo{64xlRe{cf)cf + ^(3 - 4xo)xl\cf\') 



■nil 



(19 + 121n(4^))Cfff^ + dxo{^xo{-3 - 9xo + 28xl)Re{cf)C'. 
mi. Jo o 



3Mr'eff|2 



xo(-3 + 3a;o + 48x5 - 80x^)|Q 



m 



2 ' 



(D-90) 



(-119m^2 _ i44^K)^2 _ go(^^2 _ ^12) in(4Z!^))ceff' + -^(2m| - Sm^j^i^o 



rrir 



+ 



+ 



dxo-^xl{-ml - 4m'lxo + 2m%xl + 2rnjxl)Re{Cg'''-)Cj 



ms/ruB 
l(l+m2/m|) 

■nHs/mB 



64 



45m% 
'eff^^efF 



_i(55 + 841n(44))C?ff%43^ 



TT 6 ruB o niQ 

+ dxo{^{6 - 5xo)xlRe{cf)cf + ^(18 - 39xo + 20xg)x^|G 
Jo o 9 

(11 + 41n(4^))Cfff' + 13^ + dxo{^xl{-3 - 6xo + 35xl)Re{cf)cf 
rrin yu Jo 3 



■iefF|2j 



Ai 



mr 



^x^(-3 + 6x0 + 51x^ - 100x^)|C|^|2 



m 



2 ' 



(D-91) 



(xo(so-m^)) — 
Bo 



+ 



TT 



+ (ixo(^(3 - 5xo)xlRe{cf)cf + ^(9 - 27xo + 20xg)xg|Ci 
Jo 3 9 



2 



m 



m 



-l{5 + 3ln{4^))cf +13^ 



.C 



10 



m 



90 



(D-92) 



+ r dxo(^x3(3 + 7xo)Re{cf)cf + ^x^S + 3xo - 20x2)|C; 
Jo 3 3 



|ff|2 



m 



2 ' 



(so - m^) 



'0 



= ^^(1,0)^2 + 



+ 



+ 



+ 



TT 



m 



(23 + 241n(4^))C|ff +13 



m 



B 
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30 



/ ' dxo(128(l - xo)xlRe{cf)cf + ^(3 - 7xo + 4x2)x2|C; 
Jo 3 



-^(31 + 241n(4^))cf' 

O III D 
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/ dxo{ 
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Ai 



rrir 



(D-93) 



x2(3 + lxo)Re{cf)cf + ^xg(3 + 3xo - 2Qxl)\cf\^ 
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-|-(119 + 601n(45))Ceff^-16^ 
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(a) (b) 




E„ [GeV] S„ [GeV2] 

Figure 10: Hadron spectra in B ^ Xgi'^i'in the Fermi motion model with the cuts on the dilpeton 
mass defined in eq. ^7^); (a),(c),(e) for the hadronic energy and (b),(d),(f) for the hadronic invariant 
mass corresponding to cut A,B,C, respectively. The solid, dotted, dashed curves correspond to the 
parameters (Ai,A) = (-0.3, 0.5), (-0.1, 0.4), (-0.15, 0.35) in (GeV^ , GeV), respectively. 
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